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We investigate the skyrmion formation process in nano-structured FeGe Hall-bar devices by mea¬ 
surements of topological Hall effect, which extracts the winding number of a spin texture as an 
emergent magnetic field. Step-wise profiles of topological Hall resistivity are observed in the course 
of varying the applied magnetic field, which arise from instantaneous changes in the magnetic nano¬ 
structure such as creation, annihilation, and jittering motion of skyrmions. The discrete changes in 
topological Hall resistivity demonstrate the quantized nature of emergent magnetic flux inherent in 
each skyrmion, which had been indistinguishable in many-skyrmion systems on a macroscopic scale. 

PACS numbers: 72.25.Ba, 75.25.-j, 85.75.-d 


Recent studies on condensed matters with peculiar 
band topology M have established that emergent mag¬ 
netic field arises from the Berry curvature in the mo¬ 
mentum space and leads to unique electromagnetic re¬ 
sponses |5 HlZ|, such as quantized anomalous Hall effects 
in topological insulators [MS]. Even for originally non- 
topological band (momentum-space) structures, materi¬ 
als may acquire a topological nature when the electrons 
interact with a specific spin texture in real space. The 
magnetic skyrmion in chiral magnets is an intriguing ex¬ 
ample of such topological spin textures, in which con¬ 
stituent spins point in all directions wrapping a sphere 
[Til 121. Such a winding spin texture whose topology 
can be labeled as the integer winding number —1 exerts 
a quantized emergent magnetic flux of — cf>o = —h/e on 
conduction electrons in case of strong limit of spin-charge 
coupling; as a result, topological Hall effect (THE), 
besides normal and anomalous Hall effects (NHE and 
AHE), can be observed as a hallmark of the skyrmion 
formation 13-1$). 


In addition to the topological aspect, the skyrmion pos¬ 
sesses the small-sized particle nature (its typical diame¬ 
ter is 3-200 nm), which potentially enables magnetic non¬ 
volatile memories with ultrahigh density |l9U20l]. To fully 
benefit from such nano-sized characteristics, a methodol¬ 
ogy for identifying individual skyrmion needs to be devel¬ 
oped. Although it has been demonstrated that observa¬ 
tion and even manipulation of a single skyrmion are feasi¬ 
ble by means of Lorentz transmission electron microscopy 
121 and spin-polarized scanning tunneling microscopy 
21 . [22! ) . identifying individual skyrmions by electrical 
means remains as a challenge for developing skyrmion- 
based magnetic memories. In this context, we can envis¬ 
age that the quantized emergent magnetic flux inherent 
in the skyrmion may play a crucial role in identifying a 
single skyrmion. Most notably, smaller skyrmions for in¬ 


creasing recording density facilitate their detection with 
use of the emergent flux, contrary to conventional case 
where smaller magnetic bits are generally more difficult 
to detect: Effective magnetic field observed in Hall volt¬ 
age increases inversely with the square of skyrmion size 
a s k, following the relation of _B e ff = — 4>o/a1^, which ide¬ 
ally reaches approximately 4000 T in a 1-nm skyrmion. 
Nevertheless, such a remarkable feature of skyrmions re¬ 
mains elusive, presumably because a large number of 
skyrmions involved in the macroscopic system studied 
so far inevitably smears out the quantized nature of the 
skyrmion 2^, 24 ). 


In this Rapid Communication, we demonstrate dis¬ 
cretized changes in the emergent magnetic flux arising 
from a finite number of skyrmions, by measuring topo¬ 
logical Hall resistivity of nanoscale FeGe Hall-bar devices. 
By measuring the inclination-angle dependence of topo¬ 
logical Hall resistivity, which works as a test for skyrmion 
formation in thin films [25j |. we can ensure a hysteretic 
formation of skyrmions appearing as a hysteresis loop in 
Hall resistivity [261] . In the devices with circuit line width 
between 50-250 nm, the discretized profiles of topological 
Hall resistivities show up in the hysteresis loops of Hall 
resistivities where rapid conversions between skyrmions 
with opposite core-magnetizations occur. Further minia¬ 
turization to 32-nm-wide circuit removes the typical sig¬ 
nature of skyrmion formation (a hysteretic topological 
Hall resitivity), which indicates that skyrmions are not 
able to form in a smaller area than half the size of them¬ 
selves. 


Hall-bar devices of various circuit-line-widths (32 nm 
- 250 nm and 10 fj, m) were fabricated from 40-nm-thick 
FeGe epitaxial thin film by using electron-beam lithog¬ 
raphy method. A top-view image of a 50-nm-wide de¬ 
vice is examplified in Fig. 1(a). The electrical leads for 
current and voltage were made to have the same width 
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(re); skyrmion formations in the intersectional areas of 
w x w were probed by measurements of the topological 
Hall effect [27] - The FeGe film was grown on a highly- 
resistive Si (111) substrate (p > 1000 H cm) by codepo¬ 
sition of Fe and Ge at a substrate temperature of 325 
°C by molecular beam epitaxy. A very low concentra¬ 
tion of impurity phases was verified by high-angle X-ray 
diffraction (XRD) as detailed in Supplemental Material 
[ 27 ] ]. Here we note that the electrical conduction at 2 
K is not the ballistic one because the residual resistivity 
ranges between 50-150 pH cm ( 27 }- We do not have to 
take mesoscopic effect into account. 

We first confirm skyrmion formation in our FeGe thin 
film by THE. Topological Hall resistivity p yx associated 
with formation of a non-coplanar magnetic structure is 
usually defined as an additional component to the terms 
varying linearly with magnetic field H and magnetization 
M, i.e., the normal and anomalous Hall resistivity py X 
and Py X , respectively. However, those two components 
(Py X and Py X ) may show H- and M-nonlinear dependence 
especially in multi-band metals S151- Here we adopt 
another precise method for evaluating THE based on the 
H-direction-sensitive stability of skyrmions in thin films 
[25| . According to neutron scattering experiments on 
skyrmionic materials 0, the plane of skyrmion lattice, 
where the three magnetic modulation vectors Q lie, is 
strictly perpendicular to the magnetic field. Namely, en¬ 
ergy gain from interactions among the spins on the same 
plane perpendicular to H is crucial for stabilization of 
skyrmions. When H is tilted from the normal vector of 
the thin film plane n, Q s have out-of-plane components, 
which means that spins near the surfaces lose their part¬ 
ner spins to interact with [Fig. 1(b)]. The skyrmion state 
is thus easily destabilized with H tilted from n in thin 
films, and this feature should become more significant in 
thinner films compared to their skyrmion sizes. 

Figure 1(c) shows magnetic-field dependence of Hall 
resistivity at various inclination angles 6 at 2 K in a 10- 
pm-wide Hall-bar device. We observe a loop in Hall re¬ 
sistivity at 9 = 0°, which is recognized to be induced 
by residual skyrmions of metastable state dependent on 
magnetic-field history, as reported by Huang and Chien 
(26(. The hysteresis loop dramatically shrinks with a 
slight change in the inclination angle: the loop size be¬ 
comes less than the half just tilting at 9 = 4° and nearly 
vanishes above 9 = 10°. This sudden reduction of the 
Py X loop, which is highly unlikely for the magnetization- 
related anomalous Hall effect in ordinary ferromagnets, 
represents the destruction of two-dimensional skyrmion 
structures as expected. 

We extract topological Hall component in Fig. 1(d) by 
subtracting p yx at a high angle of 6 = 30° consisting of 
normal and anomalous Hall components from low-angle 
data [ 27 J. If the vertical component of H to the film plane 
( H± ) is lower than the critical field H c , where the ferro¬ 


magnetic state is induced, the difference between Hall 
resistivities at low and high (9 = 30°) angles reads 

A Pyx = Pyx(P) Pyxi ^0 ) 

= \Pyx( e ) + pU° ) + PyM - [^,(30°) + /&(30°)] 

= {R 0 B ± + R s Ad_ i_ + R [)/Aft cos 0) — {RqB± + R S M j_) 
= R()B e ff cos 6. 


Otherwise ( H± > H c ), 

Ap y x = lp™ x (9) + Py X (9)} - [p* (30°) + p$ x ( 30°)] 
= (R 0 B± + R S M S cos 9) — ( RqB± + R S M S cos 30°) 
= R s M s (cos6 — cos 30°). 


Here Rq and R s are normal and anomalous Hall coeffi¬ 
cients, and M± and M s the vertical component of mag¬ 
netization M and the saturated magnetization, respec¬ 
tively. Although the extracted A p yx inevitably includes 
anomalous Hall component R s M s (cos9 — cos 30°) at high 
fields, A p yx for H± < H c is almost equivalent to pj jx . 
The rapid shrinkage of p yx is again highlighted in Fig. 
1(d). As magnetization at the core of the skyrmion is 
anti-parallel to the field ll|, 12], skyrmions with posi¬ 
tive (negative) core-magnetization is stabilized at neg¬ 
ative (positive) magnetic field. Hereafter, we refer to 
skyrmions with positive (negative) core-magnetization as 
core-up (down) skyrmions, which produce positive (neg¬ 
ative) emergent magnetic fields. The hysteresis loop of 
Py X demonstrates that core-up or core-down skyrmions 
are stabilized even at zero magnetic field depending on 
a magnetic-field history, as discussed later. Possibility of 
other interpretations of the hysteretic Hall resistivity is 
discussed in the Supplemental Material [2?]: nevertheless, 
the skyrmion formation is the most plausible scenario to 
explain the hysteretic behavior. 

To pursue the quantized nature of emergent flux, we 
fabricated Hall-bar micro-devices and found that such a 
hysteresis loop representing skyrmion formation remain 
robust down to a 50-nm-wide device. We show the line- 
width dependence of Hall resistivity in Fig. 2. Devices 
with wider circuit lines than 50 nm show the characteris¬ 
tic hysteresis loops representing skyrmion formation [2§| ■ 
By contrast, Hall resistivity in a 32-nm-wide circuit al¬ 
most traces the magnetization curve [ 13 ] • The collapse of 
the hysteresis loop with narrowing the circuit line repre¬ 
sents that the critical size for the formation of skyrmions 
of FeGe lies between 32 nm and 50 nm [30]; this criti¬ 
cal size is to be compared with the helimagnetic period 
(» 70 nm) or the lattice constant (« 80 nm) of skyrmion 
crystal of FeGe bulk crystal. We assume that a helix-like 
magnetic structure with no topological charge is realized 
in the 32-nm device as theoretically indicated by Du et 
al. 0 . 


Taking a close look at the hysteresis loops in the nano¬ 
scale circuits, we notice that Hall resistivity shows dis¬ 
continuous changes with field variation, which indicates 
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FIG. 1. (color online), (a) Scanning electron microscope im¬ 
age of a Hall-effect measurement circuit with line width of 
w = 50 nm. (b) Schematic illustrations of the skyrmion- 
lattice state in a 40-nm-thick FeGe thin film. The magnetic 
modulation period in the FeGe thin film is assumed to be the 
same length of 70 nm as the bulk sample. Magnetic-field de¬ 
pendence of (c) Hall resistivity p yx and (d) A p yx at various 
inclination angles at 2 K in a 10-pm-wide circuit. (See text 
for the definition of A p yx .) A sudden shrinkage of hysteresis 
loops in panels (c) and (d) represents declines in skyrmion 
numbers. 




p 0 H(J) u 0 H (T) 


FIG. 2. (color online), (a)-(f) Magnetic-field dependence 
of Hall resistivity for FeGe nano-scale circuits with vari¬ 
ous widths ranging from 250 nm to 32 nm. We subtracted 
the FI-symmetric component using the equation p yx (zkH) = 

_ 1 JVy(H)-V y {-H)]t 
21 


the discrete change in emergent magnetic field. Figure 3 
shows magnified drawings of hysteresis loops (insets) and 
difference between Hall resistivities in increasing and de¬ 
creasing field processes, A Py X (main panels), in various- 
size devices; A Py X is plotted as the normalized value by 
its maximum value in the figure. Since the conventional 
Hall components, i.e., p^ x and p^ x , do not show hystere¬ 
sis behaviors (see Supplemental Material for the M-H 
curve, which shows little hysteresis behavior [27].), A Py X 
associated with the hysteresis loop represents the genuine 
contribution from the topological Hall effect. In contrast 
to the smooth variation of A Py X in the 10-^m device 
[Fig. 3(a)], A Py X in nano-circuits of w = 50-250 nm ex¬ 
hibit step-wise profiles [Figs. 3(b)-3(f)]. The noise level 
is one order of magnitude smaller than the smallest step 
height, 5 nil cm, indicating that the discretized step can 
not be electrical noise of our experimental system [23 ]. 
Constricted geometry comparable to the size of a few 
skyrmions emphasizes skyrmions’ individuality, so that 
we successfully observed the quantized nature of emer¬ 
gent magnetic field as discrete changes in topological Hall 
resistivity mainly originating from creation/annihilation 
of skyrmions in changing the applied magnetic field. 

Incidentally, we found that the step heights display 
considerable variation and are not equal to integer times 
of unit topological Hall resistivity Py X = — Ro4>o/w 2 , 
corresponding to creation/annihilation of one skyrmion. 
This indicates additional contributions to the discrete 
changes in emergent magnetic field, which we attribute to 
discontinuous motions of skyrmions via trapping by or re¬ 
leasing from impurity or defect sites in the course of field 
changes: this bears analogy to the Barkhausen effect, i.e., 
discontinuous changes in magnetization via the similar 
motion of ferromagnetic domains. The possible fluctua¬ 
tion of skyrmions’ positions, especially skyrmions putting 
their feet on and off the verge of the probed square area 
of the Hall device, may produce “halfway” discrete topo¬ 
logical Hall effect. For examples, the fluctuation of p yx 
around p^H = 0.2-0.3 T in 120-nm-wide device [Fig. 
3(d)] and dip structure of Py X around poH = 0-0.1 T 
in 100-nm-wide device [Fig. 3(e)] can be assigned to dis¬ 
continuous changes in skyrmion position. Here we note 
that the steady flow of skyrmions via the spin-transfer 
torque, which may also contribute to the variation in the 
step heights, probably does not occur in present study 
considering the enhanced critical current density in nano¬ 
structures [13, EH. 

Hall resistivity data enable us to envisage the skyrmion 
formation process in the nano-structured Hall bars. Here 
we take, for example, the case of w = 250 nm wide de¬ 
vice, which can accommodate approximately 9 skyrmions 
at the overlap of current and voltage leads w x w. Fig¬ 
ures 4(a) and 4(b) are the conceivable development of 
numbers of core-up and -down skyrmions in increasing 
and decreasing field processes, respectively. In increas- 
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FIG. 3. (color online), (a)-(f) Magnetic-field dependence 
of difference between Hall resistiveities in increasing and de¬ 
creasing field processes, which is denoted by A p^ x . Normal¬ 
ized A p^ x by its maximum value is presented. Insets are 
magnified images of each p yx at low fields. 


ing magnetic field from a large negative field (< —H c ) 
[Fig. 4(a)], the magnetic state is anticipated to trace the 
successive stages as illustrated in Fig. 4(e): first, core- 
up skyrmions gradually pop up [panel (ii) of Fig. 4(e)] 
in the background of ferromagnetic state of pointing- 
down moments (i); next, high-density core-up skyrmions, 
possibly skyrmion lattice state, is realized and remains 
around zero field as the metastable state (iii); as the mag¬ 
netic field increase from zero to a positive value, core-up 
skyrmions are destabilized to rapidly disappear coalesc¬ 
ing with neighboring ones (iv); and then form a convo¬ 
luted helical structure (v); the helical structure fragments 
to form core-down skyrmions (vi), which then crystallize 
(vii) and subsequently begin to disappear (viii); finally, 
spins are fully polarized upward to form the ferromag¬ 
netic state (ix). 

Magnetic-field profile of Hall resistivity can be simu¬ 
lated, as shown in Figs. 4(c) and 4(d), from the assumed 
skyrmion formation [Figs. 4(a) and 4(b)]. A unit of 
topological Hall resistivity produced by a single skyrmion 
is determined from A p^ x : Since the maximum value of 
A py X represent the difference between topological Hall 
resistivities produced by close-packed core-up and -down 
skyrmions, the unit of pJ JX is A Py X divided by twice the 
maximum capacity (i.e., 2x9 skyrmions) of the overlap 
area (w x w). In the present example, the unit of Py X is 
derived as approximately 5 nfl cm from Apy X /(2x9). We 
show the simulated Py X in Fig. 4(d) along with NHE and 
AHE evaluated from the slope of Hall resistivity at high 


Conceptual diagrams 






FIG. 4. (color online). Conceptual diagrams for skyrmion 
formation and corresponding Hall resistivities in a 250-nm- 
wide nano-scale circuit. Variations of numbers of core-up and 
core-down skyrmions as functions of magnetic field in (a) in¬ 
creasing and (b) decreasing field processes. Corresponding (c) 
Hall resistivity and (d) its decomposition into p^ x (NHE), p^ x 
(AHE), and p^ x (THE), (e) Series of schematic illustrations 
representing an expected development of magnetic structure 
as a function of magnetic field. 


fields above H c and magnetization M, respectively. The 
overall profile of pj jx bears a remarkable resemblance to 
the extracted py X of the 10-^mi device [Fig. 1(d)], which 
supports our model. The total Hall resistivity as simu¬ 
lated [Fig. 4(c)] also shows a good agreement with the 
experimental results [Fig. 2(a)]. The discretized feature 
of Py X especially shows up in parts of the hysteresis loop 
where the rapid switch between core-up and core-down 
skyrmions occurs, which we focus on in Fig. 3. 

In conclusion, we observed the discretized topologi¬ 
cal Hall effect in nano-structured FeGe Hall-bar devices 
with their sizes of one or a few skyrmions. The dis¬ 
continuous changes in topological Hall resistivity, arising 
from discontinuous motion, creation and annihilation of 
skyrmions, represent the quantized nature of the emer¬ 
gent magnetic field inherent in each skyrmion, which gets 
buried in many other skyrmions residing in a large-size 
sample used in prior experimental studies. The present 
study shows that the emergent magnetic field produc¬ 
ing the topological Hall effect can be utilized to detect 
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single skyrmion in future skyrmion-based nanodevices. 
Skyrmion formations in FeGe films and nano-devices re¬ 
quire further validation by real-space observations, which 
are left as future challenges. 
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